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Abstract 

We analyze electroproduction of light vector mesons {V = p,4> and uj) at small 
Bjorken-x in the handbag approach in which the process factorizes into general parton 
distributions and partonic subprocesses. The latter are calculated in the modified 
perturbative approach where the transverse momenta of the quark and antiquark 
forming the vector meson are retained and Sudakov suppressions are taken into ac- 
count. Modeling the generalized parton distributions through double distributions 
and using simple Gaussian wavefunctions for the vector mesons, we compute the lon- 
gitudinal cross sections at large photon virtualities. The results are in fair agreement 
with the findings of recent experiments performed at HERA and HERMES. 



1 Introduction 

Recently we analyzed light vector-meson electroproduction in the generalized Bjorken 
regime [1]. That study bases on QCD factorization [2], Sj of the process 7*p Vp into 
hard parton-Ievel subprocesses - meson electroproduction off partons - and soft proton 
matrix elements representing generalized parton distributions (GPDs). The subprocesses 
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themselves factorize into hard quark-antiquark pair production off partons - amenable to 
perturbation theory - and soft qq transitions to the vector mesons. It has been shown 
[2], [3] that in this so-called handbag approach meson electroproduction is dominated by 
transitions from longitudinally polarized virtual photons to vector mesons polarized alike. 
Other transitions are suppressed by inverse powers of the virtuality of the photon, Q^. In 
Ref. [1] we examined the kinematical region accessible to the HERA experiments which is 
characterized by high energies and very small values of Bjorken's variable, xbj (<10~^). 
In this region vector-meson electroproduction is under control of the gluonic GPDs and 
the associated gluonic subprocess ^9 ~^ the quark GPDs play only a minor role. In 
the present work we are going to extend our previous analysis [1] to lower energies and 
to values of xbj up to about 0.2. This extension necessitates the inclusion of sea and va- 
lence quark GPDs into the analysis as well as the associated subprocess 7*g — > Vq. As 
it turns out from our analysis, even in the kinematical region accessible to the HERMES 
experiment, the gluonic GPD provides substantial contributions to vector- meson electro- 
production. This observation is in conflict with results from a previous attempt [1] where 
only the quark contributions have been calculated within the handbag approach while the 
gluonic one has been estimated from the leading-log approximation [5] (where the gluon 
GPD is approximated by the usual gluon distribution) and added to the quark contribution 
incoherently. On the other hand, in the recent leading-twist handbag analysis of meson 
electroproduction performed by Diehl et al [6j a relative strength of gluon and quark con- 
tributions has been found that is very similar to our result. It is to be stressed however 
that the handbag approach to leading-twist order grossly overestimates the longitudinal 
cross section in the kinematical region accessible to current experiments. 

Here in this study we will restrict ourselves to the analysis of the longitudinal cross 
section, the most important and least model-dependent observable of vector meson elec- 
troproduction. In Sect.[2]we will briefly recapitulate the handbag approach. In some detail 
we will only discuss the formation of the vector meson from a quark-antiquark pair. In 
order to cure the mentioned deficiencies of the leading-twist mechanism we will employ 
the so-called modified perturbative approach [7] in which the quark transverse degrees 
of freedom are retained and Sudakov suppressions are taken into account. In Sect. [3] we 
construct the GPDs required in the factorization formula, from a reggeized ansatz for the 
double distribution [8]. How we fix the parameters that specify the GPDs will be described 
in Sect, m where also numerical results for the model GPDs are presented. The comparison 
of our results for the longitudinal cross sections with experiment is left to Sect. [51 Our 
summary is presented in Sect. El 

2 The handbag amphtude 

We are interested in the process 7*p — > Vp for longitudinally polarized photons (72) 
and vector mesons (V^). This process can be extracted from electroproduction of vector 
mesons by exploiting the familiar one-photon exchange approximation. We work in a 
photon-proton center of mass system (c.m.s.), in a kinematical situation where the c.m.s. 
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energy, W, as well as the virtuality of the photon are large while Bjorken's variable 

is small (xbj ^0.2). The masses of the nucleon and the meson are denoted by m and mv, 
respectively. Mandelstam t is assumed to be much smaller than Q^. The proton has a rich 
structure. For the dominant parton helicity non-flip configurations in the subprocess there 
are four GPDsl for each type of partons, named H, H, E and E. As discussed in detail 
in Ref. [I] for unpolarized protons and small xbj, respective small skewness 

+ (2) 

only the GPD H is to take into consideration, the three other ones do not contribute {H, E) 
or can be neglected (E). Two subprocesses contribute to meson electroproduction, namely 
1*9 ~^ Vg and 7*g Vq (see Fig. [1]). This prompts us to decompose the •yip — > VlP 
amplitude accordingly 

Mv = M'y + M'l,. (3) 

The amplitude is normalized such that the partial cross section for longitudinally polarized 
photons reads (A is the usual Mandelstam function) 

dai 1 



dt 



\Mv\. (4) 



Strictly speaking this cross section also receives contributions from the amplitudes for 
transitions from longitudinally polarized photons to transversally polarized vector mesons. 
However, as the analysis of the spin density matrix elements of the vector mesons reveal, 
see for instance Refs. [11], |13] or [1] , this amplitude is very small and neglected by us. 
The gluonic contribution to the amplitude reads 

M'y = e^e.C^ f' dxnUx,^,Q^t' = 0)H^{x,^,t'), (5) 

while the quark one is 

Ml = eJ2eaC^ /' dxnUx,^,Q^t' = 0)H''{x,^,t'). (6) 

The sum runs over the quark flavors a and denotes the quark charges in units of the 
positron charge e. The non-zero flavor weight factors, Cy, read 

c; = -c,'^ = C = = 1/v^, c; = i. (7) 



Parton helicity flip configurations provide four more GPDs [9]. Their neglect is vindicated by the 
properties of the subprocess amplitudes which provide factors of either —t/Q^ for the gluonic subprocess or 
y/^-i/Q for the quark one p^lllOj. The latter process is further suppressed by a twist-3 meson wavefunction. 

^ In Eq. ^ also terms oc x^jm^/Q^ and oc XB^t/Q^ occur which can safely be neglected in the kine- 
matical region of interest. 
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Figure 1: Typical lowest order Feynman graphs for the two subprocesses of meson electro- 
production. 

Only the t dependence of the GPDs is taken into account in the amplitudes ([5]) and ([6]). 
That of the subprocess amplitudes 7i provide corrections of order t/Q"^ which we neglect 
throughout this paper. In contrast to the subprocess amplitudes the t dependence of the 
GPDs is scaled by a soft parameter, actually by the slope of the diffraction peak. The 
full subprocess amplitude, given for instance in Ref. fTO], indicate a breakdown of coUinear 
factorization to order t/Q"^. As the argument of the GPDs we use t' which is defined as 

t t imin 5 (8) 

where 

t^in = -4m2 , (9) 

1 - r 

is the minimal value of t allowed in the process of interest. This way we take into account 
a kinematical power correction. Since tmin is small, even at ^ ~ 0.1 it only amounts to 
-0.036 GeV"^ this power correction, absorbed into the GPD, is tiny. Also other power 
corrections of kinematical origin, as for instance given in Eq. ([2]) or in the phase space 
factor (jlj), are taken into account by us. With the exception of these kinematical effects 
hadron masses are otherwise neglected. 

Let us now turn to the discussion of the subprocess amplitudes. As is well-known, for 
the kinematics accessible to current experiments the leading-twist contribution in which the 
partons are emitted and reabsorbed by the protons collinearly and the meson is generated 
by one-gluon exchange in coUinear approximation, does not suffice, see e.g. [H E]. The 
longitudinal cross section ctl, i.e. the integrated differential cross section (jlj), calculated to 
leading-twist order is well above experiment although with the tendency of approaching 
experiment with increasing Q^. In Sect. |5]we will return to this issue and present more 
details on it. Another indication of the failure of the leading- twist order is the smallness 
of R, the famous ratio of the longitudinal and transversal cross sections. 

As is well-known from extensive studies of electromagnetic form factors at large mo- 
mentum transfer leading-twist calculations are instable in the end-point regions since con- 
tributions from large transverse separations, b, of quark and antiquark forming the meson 
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Figure 2: Illustration of gluonic radiative corrections in axial gauge. 

are not sufficiently suppressed. In oder to eliminate that defect the so-called modified per- 
turbative approach has been invented [7] in which the quark transverse degrees of freedom 
are retained and the accompanying gluon radiation is taken into account. Thus, for the 
quarks and antiquarks entering the meson one allows for quark transverse momenta, k^, 
with respect to the meson's momentum. In addition, as suggested in Ref. |14j, one also 
makes allowance for a meson light-cone wavefunction \E'v'(r, k±) where r is the fraction of 
the light-cone plus component of the meson's momentum the quark carries; the antiquark 
carries the fraction f = 1 — r. In Ref. [7j the gluon radiation has been calculated in 
the next-to-leading-log (NLL) approximation using resummation techniques and having 
recourse to the renormalization group. The quark- antiquark separation b in configuration 
space acts as an infrared cut-off parameter. Radiative gluons with wave lengths between 
the infrared cut-off and an upper limit (related to the hard scale Q"^) yield suppression, 
softer gluons are part of the meson wavefunction while harder ones are an explicit part 
of the subprocess amplitude. Congruously, the factorization scale is given by the quark- 
antiquark separation, fip = in the modified perturbative approach. In axial gauge 
the Sudakov factor can be regarded as a modification of the meson's wavefunction [7] in a 
fashion that is depicted in Fig. [2l 

Here in this work we are going to employ the modified perturbative approach too. In 
contrast to Ref. [1| we still consider the partons entering the subprocess as being emitted 
and reabsorbed by the proton collinearly. This supposition relies on the fact that all Fock 
states of the proton contribute to the GPDs. Hence, the r.m.s. of the partons inside 
the proton refiects the charge radius of the proton (i.e. (k^^)^/^ ~ 200 MeV). Only a mild 
k± dependence of the GPDs is therefore to be expected. This is to be contrasted with 
the situation for the meson where the hard process only feeds its valence Fock state. The 
compactness of the latter entails much larger values of k^. The modified perturbative 
approach applied to the subprocess, is to some extent similar to the mechanism proposed 
in Ref. [15j for the suppression of the gluon contribution to meson electroproduction in the 
leading ln(l/xBj) approximation [5]. 

Since the resummation of the logs involved in the Sudakov factor can only efficiently be 
performed in the impact parameter space [7] we have to Fourier transform the lowest-order 
subprocess amplitudes to that space and to multiply them there with the Sudakov factor. 
This leads to {i = g,q) 

= j d^d% (r, -b) Pl.{x, e, r, Q\ b) a,{^iR) exp[-5(r, b, Q^)] . (10) 

The two-dimensional Fourier transformation between the canonical conjugated b and 
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spaces is defined by 



/(b) = rf2k^expHk^-b]/(k^). (11) 

The Sudakov exponent S in ffTOl) is given by [7] 

S{t, b, Q) = s(r, 6, Q) + .(f , 6, Q) - ^ In (/^/^Mq^^) ^ (^2) 

where a Sudakov function s occurs for each quark fine entering the meson and the abbre- 
viation 

b = -ln(6AQCD), (13) 

is used. The last term in (|T2l) arises from the apphcation of the renormahzation group 
equation {(3q = 11 — |?t-/) where nj is the number of active fiavors. A value of 220 MeV for 
^QCD is used here and in the evaluation of as from the one- loop expression. The renor- 
mahzation scale is taken to be the largest mass scale appearing in the hard scattering 
amplitude, i.e. = max{rQ,fQ, 1/b). This choice avoids large logs from higher orders 
pQCD. Since the bulk of the handbag contribution to the amplitudes is accumulated in 
regions where ^r is smaller than 3 GeV we have to deal with three active fiavors, i.e. we 
take Uf = 3. For small b there is no suppression from the Sudakov factor; as b increases 
the Sudakov factor decreases, reaching zero at 6 = 1/Aqcd- For even larger b the Sudakov 
is set to zero 0. The Sudakov function s reads 



s(r, b,Q) = — l^q In + ^^L - terms , (14) 

where 

q = In (rQ/(v^AQCD)) . (15) 

Actually we do not use the version of the NLL terms quoted in Ref. [7] but rather that 
one given in Ref. [16]. The latter one includes some minor corrections which are hardly 
relevant numerically. Due to the properties of the Sudakov factor any contribution to 
the amplitudes is damped asymptotically, i.e. for In^Q"^ / Aqqj^) — > oo, except those from 
configurations with small quark-antiquark separations. 

The hard scattering kernels J-'y or their Fourier transform J-'y occurring in Eq. (fTOj) 
are computed from the pertinent Feynman graphs, see Fig. [H The result for the gluon 
subprocess is discussed in some detail in Ref. pj| and we refrain from repeating the lengthy 
expressions here. For quarks, on the other hand, the hard scattering kernel reads 




1 1 



kl + f(x + 0QV(20 - hi- r(x - 0Q7(20 - 



(16) 



^The definition of the Sudakov factor is completed by the following rules [7]: exp [—5*] = 1 if cxp [—5] > 
1, exp \-S\ = if 6 > 1/Aqcd and s(/3, 5, Q) = if 6 < \/2//3Q. 
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where Nc denotes the number of colors and Cp = {N^ ~ ^)/ (2A^c) is the usual color factor. 
Here and in the gluon kernel as well we only retain k± in the denominators of the parton 
propagators where it plays a crucial role. Its square competes with terms oc r(f)(5^ which 
become small in the end-point regions where either r or f tends to zero. 

The denominators of the parton propagators in (JT6i) and in the gluonic kernel are either 
of the type 

Ti = ,\, (17) 



or 



(18) 



' kl - d2{x ± OQ^ - le ' 
where di > 0. The Fourier transforms of these propagator terms can readily be obtained: 



T2 



+ 



ZTT 

4 ° 



/rfi bQ) 



X 



(19) 



where Kq and Hq^^ are the zeroth order modified Bessel function of second kind and Hankel 
function, respectively. 

For the purpose of comparison we also quote the leading-twist results (i.e. the limit 
k± — >• 0) for the subprocess amplitudes 



'^V{U.) 



'^V{l.t.) 



4:7ias 



fy{l/T)v[{x-0{x-^ + ie) 
1 



-1 



fvilMvC, 



+ 



X — C, + X + C, — ie 



(20) 



Here the fy denote the decay constants of the vector mesons for which we adopt the values 
fp = 209 MeV, ^ = 221 MeV and = 187 MeV from Ref. [T7]. The 1/r moment of 
the meson's distribution amplitude which represents its wavefunction integrated over k± 
(up to the factorization scale), is denoted by (l/r)y. In deriving Eq. (!20|) use is made of 
the fact that the distribution amplitudes for the vector mesons we are interested in are 
symmetric under the exchange r 



r. 



3 Modeling the GPDs 

As in Ref. pLj we construct the GPDs with the help of the double distributions invented by 
Radyushkin [8]. The main advantage of this construction is the warranted polynomiality 
of the resulting GPDs and the correct forward limit — > 0. It is well-known that at 
low X the parton distribution functions (PDFs) behave as powers 6i of x. These powers 
are assumed to be generated by Regge poles [18] and for sea and valence quarks they 
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are identified with the ususal Regge intercepts ai{0). As a consequence of the familiar 
definition of the gluon GPD which reduces to xg{x) in the forward hmit, the power 6g is 
shifted by -1 as againts the intercept of the corresponding Regge trajectory. We generahze 
this behavior of the PDFs by assuming that the t' dependence of the GPDs is also under 
control of Regge behavior (for a similar ansatz see Ref . [I9] ) . In order to simplify matters 
we consider only two trajectories: a gluon one, ag{t'), which is understood as the hard 
physics partner of the famous Pomeron trajectory playing a prominent role in soft high 
energy diffractive scattering, and a Regge trajectory for the valence quarks, a^s,i{t'). This 
trajectory represents the family of the leading Regge poles that couple to the valence quarks 
of the proton (p, u, 02 and /2 exchange in soft scattering). Since the sea quarks mix with 
the gluons under evolution it seems plausible to use the gluon trajectory for the sea quark 
GPDs as well. This assumption is supported to some extent by the observation that the 
ratios of the various sea quark PDFs over the gluon one are approximately independent of 
at low X and for ^4 GeV^. The trajectories are assumed to be linear functions of t' 
in the small-t' range 

tti = ai{0) + a[t' , i = g,Yal, (21) 

and are accompanied by Regge residues assumed to have an exponential t' dependence 
with slope parameters bi. 

Following Radyushkin jHj but generalizing to non-zero t' , we employ the following ansatz 
for the double distributions {i = g, sea, val) 

f(6 a t') - eiMMymm' j^.(g) ^i'2ni + 2) [(1 - |/3|)^ - a^]"' 

MP,a,t)-e ""^^^^ 22-^+1 r2(n, + l) (1 - |/3|)2-«+i ' ^ ^ 

where 

hgiP) = \P\gm) rig =2, 

hteM = gsea(|/?|) Sign(/?) n,ea = 2 , 

hUP) = qUf3)m) n.ai = l. (23) 

For the decomposition of the double distribution into valence and sea contribution we 
follow the procedure proposed in Ref. [9] (for earlier discussions on this decomposition see 
H [20]) and write 

fUP,cy,t') = np,a,t')e{P) ~ n-P,a,t')e{-P). (24) 

In the forward limit (^, t' — 0) this decomposition is conform to the usual definition of sea 
and valence quark PDFs. 

According to Ref. [8] the GPDs are related to the double distributions by the integral 



/•I /-i-l/^l 

Hi{x,^,t') = d(3 da6{P + ^a-x)fi{P,a,t'). (25) 

J-l J~1+\I3\ 
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For convenience we employ an expansion of the PBFsE 

j=0 



(26) 



which is particularly useful at low j3 and allows to perform the integral ( l25ll term by term 
analytically. This results in a corresponding expansion of the GPDs 



(27) 



j=0 



For gluons and sea quarks the individual terms read (e^ = 1, esea = —1, f^gj = 3 + j/2 
6g - a'gt, mseaj = mgj - 1) [8] 



15 T{mij - 2) 
W ^{mj + 3) 



X 



{[K- + m' - xf - {ml - i)(i - e){x' - e) 

} x>i 



+ 3mjj ^(1 — x) (^^ — x) 



15 T{mij-2) 
2^ T{mij + 3) 



{x?^ [(m^^. + 2){e - xf + 3m,, ^(e' - a;)(l - x) 



-K-i)(i-a(^^-e^ 

and for valence quarks {rn^^\j = 2 + j'/2 — (5vai — Q^vai''') 

3 r(mvaii-l) 



-x)| < X < ^ . 



f28) 



^^vali(a^,^,^') 



2^3 r(m,ai,- + 2) 

X [{e- x)(xr^'^ - x™-'^) + e(i - x)(xr™'^- + x^"^)} ^ > e 

- ^ ^^""^^'^"^^^^'^ [e'-a; + m.,i,e(l-a;)l -e<a:<e. (29) 



2^2 r(mvaii + 2) 
In Eqs. ( 129|) and fl28|) we use the short- hand expressions 

x+^ x-^ 



Xi 



X2 



Xj-^ 

1 + e 



i + r i-r 

The definition of the GPDs is completed by the relations 

H^{-x, e, t') = H^ix, t') , HU-x, t') = t') 

and 



i74i(x,e,t') = 



1 > X < 



(30) 

(31) 
(32) 



^ These terms may be interpreted as a set of daughter trajectories spaced by 1/2. 
Because ^soa > 1 there is a singularity at /3 = in the integral (^5)1 . As suggested by Radyushkin [S] 

sea 



this singularity is regularized by considering (j25p for i^l^a ^ principal value integral 
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The GPDs and their derivatives up to order Ui are continuous at x = ^. In the present 
work we do not consider the full evolution of the GPDs. We rather approximate it by 
that of the PDFs. Since for ^ <C a; and t' ~ the GPDs turn into the corresponding 
PDFs up to corrections of order evolution is therefore taken into account in that region 
approximately. The corrections to the PDFs are of different size - they are very small 
for the gluon and largest for the sea quarks. The approximative treatment of evolution 
receives further support by the fact that at low skewness the dominant contribution to 
the longitudinal cross section is provided by the imaginary part of the amplitude ([3]). The 
imaginary part is related to the GPDs at x ~ ^ (here differences between x and ^ are of 
order {k'j_)/Q'^) which read for ^ -C 1: 



(1 - 6g/5 - a'//5){l - dg/4 - a'//4){l - 5,/3 - a^t'/S) ' 

^ (i-W3-«;„t73)(i-W2-«;,itV2)' ^ ^ 

For sea quarks one has to replace bg by (5sea = 1 + 5g in the first equation. Hence, at least at 
small ^, evolution of the GPDs is approximately taken into account. The implementation 
of the full evolution is left to a forthcoming paper. From Eq. fl33l) one may also read 
off the so-called skewing effect [20l [21] 122] . i.e. the enhancement of the GPDs at x = ^ 
and at t' = over the corresponding PDFs taken at the momentum fraction 2^. For the 
gluons the skewing effect provides the difference between the leading-twist result for vector 
meson electroproduction and the leading-log approximation [5] for which ~ 0) 

is replaced by 2^g{2^). It is easy to see from Eqs. (127|) . (!28|) and (!33l) that the leading-log 
approximation is only valid at very low ^. One may also see from Eq. (!33|) that the use of 
the n = 1 model for the sea quarks would lead to implausibly large values for -ffsea(^, ^, i')- 
This is another reason why we prefer the n = 2 GPD model for the sea quarks. 

In Eq. fl25|) so-called D terms for gluons and sea quarks are ignored [23]. The D terms 
ensure the appearance of the highest powers of the skewness in the moments of the GPDs. 
According to the above discussion the D terms only contribute to the less important real 
part of the amplitude since their support is the region ^ x < ^. We take this in 
vindication of neglecting the D terms. 



4 Fixing the parameters 

In this section we are going to fix the parameters of our model and to present numerical 
results for the GPDs. For the evaluation of the latter we have to choose a set of PDFs and 
to expand them according to Eq. (!26|) . Let us begin with the gluon PDF. The data used 
in current PDF analyses, e.g. Refs. [211 [25l [261 1211 [SH [29], do not constrain the parton 
distributions well for /? ^ 10~^ at low [30]. This is evident from Fig. [3] where different 
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Figure 3: Left: Current gluon PDFs versus /? at a scale of 4 GeV^ [21 EHl [26l EH EHl [29] . 
For the CTEQ6M solution only the band of Hessian errors is shown. Its width is typical 
of all PDFs. The central solid line represents the gluon PDF used in this work. Right: 
The intercept 5g versus Q"^ aiW = 75 GeV. The fit 5g = 0.10 + 0.06 In QVQg is compared 
to the HERA data for p (solid circles pjj, triangles [12]) and electroproduction (open 
circles [T5]). 



versions of I3g{(3) are displayed □ at the scale 4 GeV^; the deviations diminish with increas- 
ing scale. The uncertainties in the gluon PDF matter to vector meson electroproduction. 
Via Eqs. (1221) . (!25|) and ([5]) they propagate to the cross section and lead to corresponding 
uncertainties there. In order to overcome this deficiency we adjust the low-/5 behavior 
of the gluon PDF in such a way that good agreement with the HERA data on p and 
electroproduction is achieved. Using Eqs. (jl]), ^ and one readily obtains from the 
imaginary part of the gluon contribution 

(JL oc ly^-^-^Q') , (34) 

at fixed and small x^-^; the real part does not affect the energy dependence [1]. This 
result allows for a determination of the intercept of the gluon trajectory. We stress that 
with 5sea = 1 + the sea quark contribution leads to the same energy dependence of the 
cross section as the gluon and is in so far included in flM|) . For the large energy available 
at HERA the contributions from the valence quarks are negligible. 

The power 5g has been extracted from the HERA data [HI [121 [13] on the electropro- 
duction cross section. The results are shown in Fig. [31 Since the ratio of the longitudinal 
and transverse cross section is only mildly energy dependent the difference between the 

^ Here and in the following we denote the argument of the PDFs by (3 in parallel to the definition (|23|) 
in order to avoid confusion. 
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Efluon 


strange 






6 


0.10 + 0.06L 


i + s, 


0.48 


0.48 


Co 


2.23 + 0.362 L 


0.123 + 0.0003 L 


1.52 + 0.248 L 


0.76 + 0.248 L 


Cl 


5.43 - 7.00 L 


-0.327- 0.004 L 


2.88 - 0.940 L 


3.11 - 1.36L 


C2 


-34.0 + 22.5 L 


0.692 - 0.068 L 


-0.095 L 


-3.99 + 1.15 L 


C3 


40.6 - 21.6 L 


-0.486 + 0.038 L 









Table 1: The parameters appearing in the expansion (126|) of the PDFs {L = In Q'^/Qq, 
Ql = 4: GeV^). The expansion provide fits to the CTEQ6M PDFs [21] in the range 
10-2 < /? < 0.5 and 4 GeV^ < < 40 GeV^ The powers 6 are kept fixed in the fits. 

energy dependence of the full cross section and of o"l is marginal. A straight-line fit to the 
HERA data is also shown in Fig. [3] and its parameters are quoted in Tab. [H Keeping this 
result on 6g subsequently fixed, the coefficients Cgj in Eq. fl26|) are fitted to the CTEQ6M 
gluon PDF [24] in the range 4 GeV^ < < AO GeV^ and 10"^ < [3 < 0.5. The values 
obtained for the coefficients are given in Tab. [1] as well. The resulting fit is also shown in 
Fig. [31 In the quoted range of and for f] the fit agrees very well with the CTEQ6M 
solution; it is always well inside the band of Hessian errors, see Fig. [31 Larger values of P 
are irrelevant to us since the contribution of the GPDs from the region 0.5 ^ /? to the real 
part of the amplitude is less than 0.5%. From the gluon PDF just described we evaluate 
the gluon GPD with the help of Eqs. ( |27j) and ( |28ll . For a set of skewness values it is shown 
in Fig. [H In the context of the error assessment to be executed below, we will comment 
on the implications of the other PDF solutions. 

A power Sg rising with Q^, is untypical for the Regge-pole model. Even more important, 
Sg > leads to a cross section that increases as a power of the energy, see Eq. flMl) . At very 
high energies, it will therefore violate the Froissart bound and, hence, unitarity. Obviously, 
there must occur a saturation at some scale of that will limit the rise of the gluon PDF 
and GPD and will restore unitarity. In other words, the description of diffraction by a 
Pomeron-type pole with an intercept, ag(0), larger than unity is to be considered as an 
effective parameterization that holds in a finite although possibly large range of energy 
|31j . Unitarity will ultimately force the generation of a series of shielding cuts [32] that 
will prevent the violation of the Froissart bound. An effective parameterization may have 
parameters that dependent on the process and on the kinematics. 

For the slope of the gluon trajectory we take the value a'g = 0.15 GeV^^ which is 
slightly smaller than that of the usual soft Pomeron [31] but agrees with the value observed 
in photoproduction of the J /"if [33] and other vector mesons [31]. Small values of the gluon 
and soft Pomeron slopes are required since the diffraction peaks show little shrinkage. 

Using Ssea. = I + Sg we perform an analogous fit of the strange quark CTEQ6M distri- 
bution (assuming s(/?) = s{j3)). The parameters are quoted in Tab.[Tltoo and the resulting 
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GPD if^gg^ constructed through the double distribution fl22|) is shown in Fig. HI For the 
CTEQ6M PDFs [23] the u and d distributions at low (3 are very close to each other and en- 
hanced by an approximately Q^-dependent but /5-independent factor Ks{Q^) as compared 
to the strange quark PDF. In an attempt to keep the GPD model simple we therefore 
assume 

-^sea — ^sea — -f^sea 5 (35) 

where the flavor symmetry breaking factor is parameterized as 

Ks = 1 + 0.68/(1 + 0.52 \nQ^/Ql) , (36) 

as obtained from a fit to the CTEQ6M PDFs. Eq. fl35|) is a simplification which as one 
may object, is unjustified given the high level of accuracy the current PDF solution have 
reached. However we are constructing model GPDs which implies a theoretical uncertainty 
of unknown strength. It seems premature in the present state of the art to transfer the full 
complexity of the current PDFs to the model GPDs. It is not probed by the present data 
on vector meson electroproduction |^ and would rather confuse than elucidate the physical 
interpretation. 

For the valence trajectory we adopt standard soft physics parameters, avai(O) = 0.48 
and a'^^i = 0.9 GeV~^. This is in fair agreement with the low-x behavior of the valence 
quark PDFs at low factorization scale (up to at least 40 GeV^) within errors. Only a very 
mild, negligible effect of evolution is to be observed for these parameters. Keeping again 
the intercept fixed we fit the expansion fl26l) to the u and d valence quark distributions 
of the CTEQ6M solution and evaluate the correponding GPDs from Eqs. fl27j) and fl29l) . 
The obtained parameters are quoted in Tab. [T] and if^^j at t' = is displayed in Fig. HI 
In contrast to the sea quark GPDs where the n = 1 model is to be rejected because of its 

® Since as yet there are only data on p production available but not on lo only the combination H^^^ — 
Cd H^^^ is probed. 
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very strong skewing effect (and because of the Pomeron-type interpretation), the skewing 
effect for the valence quarks is much weaker. 

Now having specified the parameterization of the Regge trajectories we turn to their 
determination of the residues, i.e. of the slope parameters 6j. As we repeatedly mentioned 
vector meson electroproduction behaves diffractively, i.e. its differential cross section de- 
crease exponentially with t' . In the handbag approach this behavior is to be incorporated 
in the gluon GPD. Although its t' dependence appears to be more complicated than an 
exponential as is for instance seen from Eq. fl28l) or fl33|) . this is not the case in reality. They 
actually behave as exponentials as can be seen from Fig. [5] where we display H^{^,^,t') 
versus f for selected values of W and Q^. As a consequence of the smallness of a'g the 
effective slope of the gluon GPD falls together with the Regge exponential to a very high 
degree of accuracy, for the sea quark GPD the situation is similar. Hence, because of the 
repeatedly mentioned dominance of the imaginary part of the gluon contribution at HERA 
kinematics the slope of the differential cross section is given by 

5v~26, + 2a;ini^, (37) 

in our model GPDs. Insertion of Eqs. ([T]) and ([2]) translates this relation into 

High-energy data for the t' dependence only exist for the unseparated cross section da = 
dcTT + £ dai where e is the ratio of longitudinal to transversal polarization of the virtual 
photon. Ignoring possible differences between the slopes of the transverse and longitudinal 
cross sections (which can only emerge from the subprocess amplitudes) we fit bg against 
the HERA data for p [II] and production [13]. We find that the experimental slope 
parameters are well described by 

2 

bg = 6sea = 2.58 GeV-2 + 0.25 GeV"' In -— , (39) 

+ 

see Fig. [5l It is to be stressed that the term oc a'g in Eq. ( l38l) is a consequence of the 
Regge behavior while the log term in Eq. (!39|) is an ansatz. The p and slopes practically 
fall together at HERA energies; there are only minor differences at low which are not 
shown in Fig. [5l 

The zero-skewness limit of the valence quark GPDs reads 

HUx, e = 0, t') = e(''-+<a> ''^(i/^))*' qU^) , (40) 

as readily follows from fl22l) and fl25|) . This is very close to an ansatz advocated for in 
Ref. in order to extract the zero-skewness GPDs from the nucleon form factor data 
(see also Ref. [2S])- Our Regge exponential appears as the small x approximation of the 
exponential exploited in Ref. [35]. In the range 0<x<0.15 the difference between both 
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Figure 5: Left: The t' dependence of the gluon and the w-valence GPDs a.t W = 
75 GeV, = 4 GeV^ (sohd hnes) and = 5 GeV = 3 GeV^ (dashed hnes). Right: 
The slope By of the differential cross section for p and (j) electroproduction versus at 
W = 75 GeV. Data are taken from Ref. [llj (solid circles) and [13] (open triangles). The 
solid line represents the fit (l38l) . fl39|) for the case of the 0. 

the exponentials is less than 10%. The comparison with the zero-skewness analysis further 
reveals that the slope of the Regge trajectory suffices to specify the f dependence of the 
valence quark GPDs. Consequently we set 6vai equal to zero. It is also checked by us 
that the valence quark GPDs fl27|) . fl29l) respect the sum rule for the Dirac form factor at 
any value of skewness; small deviations between the form factor data and the sum rules 
evaluated from our GPDs occur at larger t' and amount to about 3% at t = —0.6 GeV^. 

Effectively the valence quarks GPDs behave as exponentials in t' too, see Fig. O The 
decisive difference is however that the valence quark GPDs show strong shrinkage due to 
the large value of a'-^^^. While at HERMES kinematics the effective slopes of the gluon and 
the valence quark GPDs are similar, is the latter much larger at HERA kinematics. 

Last not least we have to specify the meson wavefunction occurring in Eq. fllOp . As in 
our previous work pQ and in other applications of the modified perturbative approach, e.g. 
[m [16], we use a Gaussian wavefunction 

^v.(r,k^) = 87r^^Jvalexp[-a'ykl/{rf)]. (41) 

Transverse momentum integration of it leads to the associated distribution amplitude which 
represents the soft hadronic matrix element entering calculations in the collinear factoriza- 
tion approach. Actually, the wavefunction (j4T|) leads to the so-called asymptotic form of 
the meson distribution amplitude 

$As = 6rf . (42) 

Its 1/r moment occurs in the leading-twist result (1201) and acquires the value 3. The 
transverse size parameter ay is considered as a free parameter fitted to the data on the 
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integrated cross sections ctl for p and production, see next section. It can be varied 
within a certain range of values determined by the requirement that the corresponding 
r.m.s. k_L being related to the transverse size parameter by 

^ ^^av]-\ (43) 

acquires a plausible value consistent with our assumption of taking into account only the 
transverse momenta of the quarks forming the meson. A possible evolution of the transverse 
size parameter is ignored. 



5 Results for the longitudinal cross sections 

The full amplitudes for vector meson electroproduction ([5]), ([6]) are a coherent superposition 
of contributions from the various quark flavors and from the gluon. In order to shed light 
on the relative importance of the various terms we quote the leading-twist result 



» , Sttq;, 
Ms = e- 



Mn = e f, 



The integrals in Eq. (14^ read 



(x + ^)(x - ^ + ie) ' 



(x + ^)(x - ^ + ie) 

The imaginary parts of these integrals are just —TTHi[^,^,t') and for ^ ^ they exhibit 
typical Regge phases as a consequence of the analytic structure of the integrals and the 
symmetries of the GPDs (!3Ti) . 

Within the modified perturbative approach the amplitudes have the same structure as 
in Eq. fH4|) . Only the integrals fHSl) are much more complex, they do not factorize into 
products of integrals over the wave functions and such over the product of GPDs and 
propagators. Despite this the suppressions induced by the modified perturbative approach 
do not change much the relative strengths of the various contributions. One may therefore 
get an quick insight into the relative strength of the various terms from Eq. (H4l) . 

Numerical evaluation of the amplitudes ([5]) and (E]) reveals that, for skewness less 
than about 0.01, the gluon and sea contributions are dominantly imaginary while, for 
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Figure 6: Real over imaginary part of the full amplitude for p production versus W at 
= 4 GeV^ and t' = 0. 



^ ~ 0.1, their real parts are nearly as large as their imaginary part. The valence quark 
contribution behaves oppositely. The energy dependence of ratio of the real and imaginary 
parts of the full p production amplitude is shown in Fig. [6] at = 4 GeV^. A remark 
concerning the values of is in order. As an inspection of our handbag amplitude reveals 
almost the entire contribution is accumulated in a comparatively narrow region of a^- For 
instance, at = 4 GeV^ and W = 5 GeV 90% of the amplitude comes from the range 
0.3 ^ «s ^ 0.5 with a mean value of about 0.4. Hence, our handbag approach is theoretically 
self-consistent in so far as contributions from soft regions where is larger than, say, 0.6 
and where perturbation theory breaks down, are strongly suppressed. 

We are now ready to present our results for vector meson electroproduction. They 
are obtained by adjusting the transverse size parameters ay appropriately. The best fits 
provide the values 

= 0.70 GeV"^ , Up = 0.75 GeV"^ . (46) 

The corresponding r.m.s. is about 500 MeV. This value is more than a factor of 2 larger 
than the corresponding value for the quarks inside the proton and is in so far consistent 
with our assumption of taking into account only the transverse momenta of the quarks 
forming the meson. In Figs. [7] and [8] we compare our results for and p production to 
the HERA data [IH [HI [131 EZl EHl ES]- In the left panels of these figures we show 
the decomposition of the cross sections into the various contributions from the gluon and 
the quarks. The gluon contribution is dominant but the corrections from the gluon-sea 
interference amount to about 25 (50)% for (p (p) production; the valence quark contribution 
(including its interference with the gluon and the sea) to the p production cross section 
is tiny and can be neglected. The larger sea quark contribution for p production follows 

If not quoted explicitely in the experimental papers the longitudinal cross section is evaluated by us 
from information given therein. Statistical and systematical errors are added in quadrature. If necessary 
data are rescaled in W using Eq. ([M)) or exploiting the handbag predictions for the dependence on it. 
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Figure 7: The longitudinal cross section for (p production aX, W = GeV. Data are taken 
from [13] (open triangles), [37J (solid squares) and [38] (open squares). Left: Full (dashed, 
dash-dotted, dotted) line represents the handbag predictions for the cross section (gluon, 
gluon-sea interference, sea contribution). Right: Predictions for the cross section with 
error bands resulting from the Hessian errors of the CTEQ parton distributions (full line) 
and compared to the leading-twist result (dashed hne). 




[GeV^] q2 [GeV^] 

Figure 8: The longitudinal cross section for p production at = 75 GeV. Data are taken 
from [m [39] (solid squares) and [12] (open squares). For other notations cf. Fig. [7l 
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Figure 9: The longitudinal cross section for (p (left) and p (right) production versus at 
W = b (sohd line) and 10 GeV (dashed line). Data at = 5 GeV, shown as solid circles, 
are taken from HERMES [H] (for 0, preliminary) and [42] (for p). The open triangle 
presents the E665 [43] data point aXW = lQ GeV. For other notations cf. Fig. [7l 

from the flavor symmetry breaking factor Kg. In the right panels of Figs. [7] and [H] we 
display in addition to our full results their uncertainties due to the Hessian errors of the 
CTEQ6 PDFs and for comparison the leading-twist result evaluated from the asymptotic 
distribution amplitude (H2!) . Results for the cross sections evaluated from sets of PDFs 
other than CTEQ6 also fall into the error bands in most cases (an exception is set for 
instance by the PDFs determined in Ref. [29]) provided these PDFs are treated in analogy 
to the CTEQ6M set, i.e. they are fitted to the expansion fl26|) by forcing them to behave 
Regge-like with powers 5i as described above, and if necessary readjusting the transverse 
size parameters. Straight-forward evaluation of the GPDs from the various sets of PDFs 
and fixed transverse-size parameters lead to cross sections which differ markedly stronger 
than the error bands indicate. For examples see Ref. [6]. The results obtained with the 
modified perturbative approach are in remarkable agreement with the HERA data while 
the leading-twist results are clearly in excess to experiment with a tendency however of 
approaching the data and the predictions from the modified perturbative approach at 

~ 40 GeV . This in turn tells us that the effect of the transverse quark degrees of 
freedom in combination with the Sudakov suppressions become small for such values of 

while being very important at lower . Pertinent observations have also been made 
by Ivanov et al. [ID]. In their next-to- leading order leading- twist calculation of vector 
meson electroproduction large perturbative logs occur which partly cancel the leading- 
order term bringing the leading-twist result closer to experiment. These logs are included 
in the Sudakov factor ([T^ . 

In Fig. we show the results for cJi at = 5 and 10 GeV and compare them to 
the data from HERMES HH |12] and the FERMILAB experiment E665 Again we 

observe good agreement with experiment. The slopes of the differential cross section are 
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Figure 10: The differential cross section for p production versus t' aX W = b GeV and 
Q2 = 4 GeV^ 

somewhat smaller at lower energies than those at the HERA energy shown in Fig. [5l For 
instance aX W = h GeV and = A GeV^ we obtain 5.0 GeV~^ for p production and 
4.8 GeV"^ for the case of the 0. As yet the HERMES collaboration has only provided 
preliminary results for these slopes: 6.32 ± 0.72 GeV~^ at = 3.7 GeV^ for p [H] and 
4.6 ±1.2 GeV~^ averaged over the range 0.7 < < 5 GeV^ for production [15]. A slope 
for p production that is considerably larger than that for production is difficult to get in 
the handbag approach. Although it seems tempting to assign such an effect to the valence 
quark contribution by choosing a non-zero value for 6^^! this is likely not the solution since 
it would lead to a Dirac form factor that drops too fast with t. An alternative possibility 
seems to change the t' dependence of the sea quark GPD by choosing a value for bsen larger 
than that for bg which is theoretically not forbidden. On account of the different weights 
of the sea contribution in both the processes of interest this may generate a somewhat 
larger slope in the case of the p. In regard to the present experimental situation we leave 
this question unanswered for the time being. Pertinent future data may settle this issue. 
In order to demonstrate how close the t' dependence of the differential cross section to an 
exponential behavior is we show in Fig.[TO]the cross section for p production at = 5 GeV 
and = 4 GeV^ as an example. 

Consulting Fig. [TT] one sees that also the W dependence of the longitudinal cross section 
is correctly described within the handbag approach. Note that o"l in Figs. O [H [9] and [12] 
(left) is obtained by integrating the differential cross section over the range of t' that is 
used in the various experiments. Thus, we integrate from to —0.6, —0.5, —0.4 GeV for 
HERA, COMPASS (and E665), and HERMES kinematics, respectively. For cxl in Figs.[II] 
and [12] (right) on the other hand we integrate up to —0.5 GeV^ throughout but compare 
with actual data. The cross sections exhibit kinks at about W = 10 GeV, rather markedly 
for production, milder in the case of the p. These kinks are related to the sharp fall off 
of the gluon and sea quark GPDs with increasing ^ (note that ^ oc l/W"^ at fixed Q^), see 
the gluon PDF shown in Fig. [3] The additional valence quark contribution in p production 
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Figure 11: The longitudinal cross section for (j) (left) and p (right) electroproduction versus 
W aX = 3.8 GeV^ and 4 GeV^, respectively. The handbag predictions are evaluated 
from the interval —t' < 0.5 GeV^. Data for (p production are taken from HERMES [H] 
(solid circle), ZEUS |T3] (open triangles) and HI [37j (sohd square). The data for p 
production are taken from HERMES [42] (solid circles), E665 [43] (open triangles), ZEUS 
[T2] (open square) and HI [TT] (solid square). The dashed (dash-dotted, dash-dot-dotted) 
line represents the gluon (gluon + sea, (gluon + sea)-valence interference plus valence 
quark) contribution. For other notations cf. Fig. [71 




Figure 12: Left: The ratio of the longitudinal cross sections for and p production. Data 
are taken from HI [HI [37] (solid squares), ZEUS [12], [13] (open squares) and HERMES 
[m US] (solid circles). The solid (dashed) line represents the handbag predictions at 
W = 75(5) GeV. Right: Predictions for cu electroproduction versus W at Q'^ = 3.5 GeV^. 
For comparison the full cross section for uj production, measured by ZEUS [19], is also 
shown (solid diamond). 
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mitigates the kink. For p production we also show in Fig. [TT] the individual contributions 
from the gluons, sea and valence quarks. At ^ 5 GeV the latter contribution (including 
the interference with the gluon and sea quarks) amounts to about 50% of the full result 
but decreases rapidly with increasing W . It only contributes about 10% at ly = 10 GeV 
and is negligible at the HERA energy. 

All the results we presented so far are evaluated from the Gaussian wavefunction pTI) . 
One may wonder what the consequences of other choices of the wavefunction are. In order 
to provide a partial answer to this question we multiply the wavefunction fHTj) for the p 
meson with the factor 

50/81 

cT{2t-1). (47) 

It corresponds to the first two terms of the expansion of the meson distribution amplitude 
upon the Gegenbauer polynomials C^''^, the eigenf unctions of the evolution kernel for 
mesons [l6]. For an estimate of the impact of that factor on the cross section for p 
production we adopt the value 

S2(/io = 1 GeV) = 0.18, (48) 

obtained from QCD sum rules for the expansion coefficient |47j. From this value of B2 
we find that the cross section for p production increases by approximately 14(24)% at 
Q2 = 4 (20) GeV^ in the entire range of energy we examine. This is well within the 
uncertainties of our approach which are represented by the error bands shown in the various 
figures. The gradual decrease of the second Gegenbauer term with increasing scale is to 
a large extent compensated by the diminishing suppression through the quark transverse 
momenta and the Sudakov factor. In a calculation to leading-twist order the effect of the 
higher Gegenbauer terms is more pronounced at low Q^. For instance, at = 4 GeV^ 
the p cross section increases by 31% to leading- twist order using fHS]) and pp = Q. 

In the limit of negligible valence quark contributions one has the following relative 
strength of the three cross sections: p:to':0 = l:9:2/9upto fiavor symmetry break- 
ing effects as the differences in the decay constants and in the wavefunction or distribution 
amplitudes or the flavor symmetry breaking factor (l36l) in the sea GPDs. While the latter 
two effects disappear for — > 00 due to evolution, the differences in the decay constants 
are scale independent. Hence, for very large scales the handbag approach predicts 

M<P)/^l{p)^1 (I) = 0.248. (49) 

An analogous result holds for the ratio of the u and p cross sections. The HERA data 
[TTl [T2| [T3I [37] shown in Fig. [121 cire not far from the symmetry limit 2/9 especially at 
larger values of but the face values are clearly below it. The evolution effect of a p 
wave function broader in r than the Gaussian given in Eq. ( HTl) is too mild if one accepts 
the QCD sum rule estimate of the Gegenbauer coefficients, see above discussion. It does 
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not explain the dependence of the data. Thus, we are compelled to conclude that the 
bulk of the effect seen in the HERA data is due to flavor symmetry breaking in the sea. 
Indeed with Ks as given in Eq. ( !36l) we obtain the results for the ratio of the and p cross 
sections shown in Fig. [12] which agree fairly well with experiment. The cross section ratio 
at the HERMES energy is smaller as at the higher HERA energy which as a glance at Eq. 
dH]) reveals, is to be assigned to the additional valence quark contribution to the p cross 
section at = 5 GeV (see also Ref. |1H])- Indeed our results are in agreement with the 
HERMES data HHIig, see Fig. [HI 

Predictions for uj production are presented in Fig. [12] as well. They are obtained by 
assuming = a p. Only the unseparated cross section has been measured by the ZEUS 
collaboration [19] at = 70 GeV and = 3.5 and therefore we cannot directly compare 
with our results. However, assuming that the ratio of the longitudinal and transversal cross 
sections is about 2 at this kinematics, one expects a longitudinal cross section that amounts 
to about 2/3 of the full cross section and this would be in agreement with our result. The 
valence quark contribution is stronger for u than for p production as is expected from Eq. 
dH]). For our kinematical point of reference, W = 5 GeV and = A GeV^, for instance, 
it amounts to about 65% of the u cross section. Data on the longitudinal cross section 
for UJ production would be highly welcome. They would provide information on a second 
combination of the u and d quark GPDs + H'^. 

6 Summary 

We have investigated light vector-meson electroproduction within the handbag approach. 
The partonic subprocesses are treated within the modified perturbative approach and the 
GPD H for gluons, sea and valence quarks is constructed from the CTEQ6M PDFs through 
double distributions using a Regge-inspired t dependence. The GPDs H respect all theo- 
retical constraints, i.e. the reduction formulas, positivity and polynomiality as well as the 
sum rule for the Dirac form factor of the proton. From our approach we have obtained 
a fair understanding of the longitudinal cross section over a large range of energy (from 
W = 5 till 100 GeV) and photon virtualities ( from 2.5 to 40 GeV^) provided a^Bj is small. 
A remarkable outcome of our investigation is that the gluon contribution play an impor- 
tant role over the entire range of energy we have explored. The sea quarks also contribute 
considerably at all energies although to a lesser extent than the gluons. The valence quarks 
on the other hand are only of importance for p and uj production and for energies less than 
about 10 GeV. 

We have simplified the parameterization of the sea quark GPDs, not all details of the 
PDFs are transfered to them. Thus, we have reduced the sea quark GPDs to a single 
function allowing for differences only through a flavor symmetry factor. At the present 
stage of our knowledge on the GPDs a more refined model would provide a pseudo accuracy 
that does not meet the uncertainties of the GPD model and that would rather confuse than 
elucidate. With our sea quark GPDs it becomes evident that the dependence of the 
(f) - p ratio of the longitudinal cross sections at HERA energies is generated by the flavor 
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symmetry breaking in the sea. Data for uj production may perhaps force us to improve the 
parameterization of the sea quark GPDs. The data for p production alone only probe the 
combination e^if^a — tdHtca.- 

In our previous work we have also calculated the amplitudes for other transitions 
from the virtual photon to the vector meson within the handbag approach. With them we 
have achieved a fair description of the transverse cross section and the spin density matrix 
elements at HERA energies where the gluon contribution dominates. The application 
of this approach to the spin density matrix elements at lower energies necessitates the 
inclusion of the quark contribution into the analysis which is left to a forthcoming paper. 

The handbag approach with the reggeized GPDs bears similarities to other theoretical 
models for vector-meson electroproduction which are mainly applied to the high energy 
and/or very low xbj regime. The gluonic subprocess, see Fig. [H with the accompanied 
proton-gluon vertex function forms the basis of many models. The various approaches 
essentially differ in the treatment of that vertex function (Pomeron residue [211 ED] or 
the gluon PDF in the BFKL color dipole model [M] and in the leading-log approximation 
[151 [21]) and in the assumptions which of the partons in the Feynman graphs are considered 
as soft, i.e. as quasi on-shell, and which as hard. In the handbag approach with the QCD 
factorization theorems [21 [3] as foundation, the protons emit and reabsorb quasi on-shell 
partons and the associated vertex function, a soft proton matrix element, is regarded as 
a GPD. The quark and antiquark entering the final state meson are also considered as 
on-shell particles. The associated soft qq ^ V transition is parameterized as a light-cone 
wavefunction or distribution amplitude. All other partons in the graphs shown in Fig. [1] 
are highly virtual. The advantage of the handbag approach is that once the GPDs are fixed 
other hard exclusive processes, as for instance deeply virtual Compton scattering, can be 
predicted. 
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